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Abstract 



I report on some work in progress on the dynamics of extended objects 

in field theories after a rapid phase transition, as is relevant in the early 

Universe. An analytic technique, originally introduced to approximate the 

dynamics of topological defects in condensed matter systems, is extended 

to cover relativistic objects. An exact formula for the area density of 

Oh, domain walls in D dimensions is derived. 

(D 

' The realisation that spontaneous symmetry breaking in field theories might 

actually be a dynamic phenomenon in the early Universe brought with it the 
^ . possibility that topological defects were formed in the initial stages of the Big 

^ I Bang P|. The study of the formation and evolution of topological defects there- 

fore has great implications for particle physics and cosmology, as well as being 
of interest in itself as a difficult problem in the non-linear dynamics of fields. 
The fundamental picture that has emerged is that networks of topological de- 
fects exhibit self-similar behaviour |9|, The network is described by a length 
scale ^{t), which can be thought of as the average curvature radius of the p- 
dimensional surface, and this scale increases as a power of time. In the context 
of the early Universe, this power seems always to be 1, with a constant of pro- 
portionality of order 1. Hence the defects straighten out as fast as causality 
will allow. 

In this talk I report of some work on an analytic approach to the dynamics 
of defects formed in a phase transition in the early Universe Q . The programme 
is as yet in its early stages, but already has produced interesting results for the 
dynamics of domain walls, explaining the results of numerical simulations in a 
much more detailed way than the qualitative arguments that went before. 

The analytic approach is based on a method which made its first appearance 
in the condensed matter literature [Q, . The defects are replaced by scalar fields 
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u {x) specially designed so their zeros track the positions of the defects. The re- 
sulting equations of motion are simpler to deal with than the original field equa- 
tions, particularly in the case of gauge field theories. The initial condition for the 
field are selected from a Gaussian distribution, with mean U^{r]i) = (^u^{rii,:x.)'^ 

and two-point correlation function C{r]i,r) = (^u'^{-K)u^{-K + r)6AB/N'^- This 
is supposed to represent an initial high-temperature thermal state. One can 
then average the equations of motion for the u fields, to derive a linearised ap- 
proximation, which is exactly soluble. Important quantities such as the defect 
density can then be found as a Gaussian average over the appropriate functions 
of u. The result is always an expression involving powers and derivatives of C, 
or its unequal time version C(?7, r/', x, x') = (n(r/, x')tt(ry', x')). 

Consider p-dimensional defects in D spatial dimensions ( "p-branes" ) . We 
define with A'' = D—p real scalar fields, which vanish precisely at the coordinates 
of the defects where fj, takes the values 0,...,D and a the values 

0,... ,p. Thus we begin with the N equations 

u^iXf") = 0. (1) 

Differentiating once with respect to the world- volume coordinates a", we find 

dpX^'d^u^iX) = 0. (2) 

Thus the N vectors dnU'^^X) are spacelike normals to the p-brane. 

Defects always sweep out surfaces of zero average extrinsic curvature, and 
thus the surfaces of constant must also obey this condition. The equation 
that the fields satisfy turns out to be 

n'^'^(5^a,-r^,a>^ = o, (3) 

where F^j^ is the affine connection in the background space-time, and H^'^ is 
the transverse projector onto the defect worldsheet. 

At this point we shall specialise to the case N = 1, domain walls, so that 
the algebra does not obscure the central point. The equations of motion may 
be rewritten in a form which aids the averaging procedure: 

This equation is linearised by replacing terms of the third degree by a Gaussian 
two-point correlation function multiplied by u. For example, 

d^'ud^'udpu {d^'ud^u) dpU + cyclic. (5) 

We shall suppose that the Universe is described by a flat Friedmann- Robertson- 
Walker space-time, for which the affine connection is Fj[^^ = {5'^8^ + 5'^5^^ — 
gn,yg'^){a/a). Here, a{rf) is the scale factor and rj is conformal time. The 
linearised We define the following correlators 

T{r]) = {u{7], x)u{r], x)) , S{7]) = ^ {Vu{rj, x)Vu{r], x)) . (6) 



{d^d.u - rp^dpu) = 0. (4) 
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and note that C(r/) = C{r],0). It can be shown that the hnearised equations of 
motion then take the form 



U+!Mu-y^^\ = 0^ (7) 

V 

where fj,{r}) and v are functions of D. For Friedmann models, one can show that 

Mr?) = -2v{S/S) + a{r^)[D-3{T/S)], (8) 
= [D-l-iT/S)]/D, (9) 

where a{r]) = rjd/a. 

We expect scahng solutions to have S,T oc rj^, and so as long as we are 
not near a transition in the equation of state of the Universe (such as that 
between the radiation- and matter-dominated eras), /i and v'^ are constant. 
Thus, imposing the boundary condition that u be regular as r/ ^ 0, (^ has the 
simple solution 

where (1 — /i)^/4 = . The form of the initial power spectrum Pi(A;) = |itk(?7i)P 
is taken to be white noise, behaving as at small fc, as appropriate to an 
uncor related field. 

We can substitute the solution back into the correlation functions C, S, 
and T. A self-consistent solution for ^ turns out only to be possible if we take 
u = {fj. — l)/2, so that r]{S/S) = —{D + 2). We may thus obtain an algebraic 
equation for with a = rj{a/a) as a parameter. For example, in Minkowski 
space {a = 0) we find /i = 2{D + 2). We can now calculate anything that 
can be expressed in terms of local functions of the field and its derivatives, by 
averaging over the Gaussian probability distribution for the fields. The easiest 
quantity to calculate is the comoving area density of the walls. 



A = / d'^a^dD+iix - X{a)). (11) 



Making the coordinate transformation from to (fT",u) near the wall, this 
can be rewritten as 

A = S{u)\du\. (12) 
The probability distribution takes the form 



P{df^u, u) = J\fex.p 

where 




(13) 



X = T/SC-f, Yf" = -^S^{l/Sj), Z°° = l/57, Z'^ = 5'^{l/S), 

and 7 = ~ i^) • "^^^ normalisation of the probability distribution is given 
hy Af = (27r)^^/^(CTS'^)~"^/^. We are now in a position to evaluate the average 



3 



area density. A long calculation shows that 



[Y r(p + 2)/2) 

^•^^ \ 7tCT{{D + 1)/2) 

Recalling the behaviour of S and C with conformal time, we immediately 
see that we reproduce the correct scaling behaviour for the area density, A oc 
rj^^. But this technique also gives the coefficient of proportionality, which we 
can compare with numerical simulations. The comparison in D = 3 is quite 
satisfying For example, simulations ^, ^ give A ~ 1.57/"^ in a radiation 
dominated Universe, while the theoretical prediction is 2.1rj~^. In D = 2 the 
agreement is also good, although there is evidence for a small deviation from 
the ri~^ behaviour in D = 2 [^, |^. However, it should be borne in mind that 
the theoretical approach involves an approximation which is not well-controlled: 
replacing correlators by their Gaussian averages. This is not an expansion in a 
small parameter, so assuming that the higher order corrections are small. 
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